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We calculate the static and dynamic properties of single crystal, single molecule magnets consisting 
of equal spin S — 1/2 or 5/2 dimers. The spins in each dimer interact with each other via the 
Heisenberg exchange interaction and with the magnetic induction B via the Zeeman interaction, 
and interdimer interactions are negligible. For antiferromagnetic couplings, the static magnetization 
and specific heat exhibit interesting low temperature T and strong B quantum effects. We calculate 
the frequency spectrum of the Fourier transform of the real part of the time autocorrelation function 
Cii{t) for arbitrary T, B, and compare our results with those obtained for classical spins. We also 
calculate the inelastic neutron magnetic dynamical structure factor S{q^,ui) at arbitrary T, B. 
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INTRODUCTION 

Recently, there has been a considerable interest in the 
physics of single molecule magnets (SMM's), or magnetic 
molecules. 0, |, 1 1, |, | 0, | g 0, 0, Q These con- 
sist of small clusters of magnetic ions embedded within 
a non-magnetic ligand group, which may crystallize into 
large, well-ordered single crystals of sufficient quality for 
neutron scattering studies. Usually, the spins within a 
single molecule interact mainly via the Heisenberg ex- 
change interaction. In the simplest SMM's, V2, Cu2, 
two examples of Yb2, Cr2, and four examples of Fe2, the 
magnetic cores of the molecules consist of dimers of spin 
S=l/2 Cu2+, or Yb3+, S* 3/2 Ct^+ , or S = 5/2 
Fe^"*" spins, respectively. |[ |[ |l^, |l^ Low-field magne- 
tization, nuclear magnetic resonance (NMR) and electron 
paramagnetic resonance (EPR) experiments, and zero- 
field inelastic neutron scattering experiments were made 
on some of these dimers. |, |, |, ||, |l^, |l|, || Theo- 
retically, zero-field results for the time dependence of the 
autocorrelation function of the quantum 5* = 1/2 and 
5 = 5/2 dimers were presented. ]T9| 

Here we study the simplest model of interacting 
Heisenberg spins in a magnetic field H. We assume only 
two spins, which interact with each other via the ordi- 
nary Heisenberg exchange interaction J, and also with a 
constant magnetic induction B induced by the applica- 
tion of H. For simplicity, we limit our discussion to the 
spin values S — 1/2 and S = 5/2. 

We first evaluate the static magnetization M and spe- 
cific heat Cv as functions of T and B. We find that for 
ferromagnetic (FM) exchange couplings J > 0, both M 
and Cv behave at low T as do a single paramagnetic ion 
with spin 25, qualitatively similar to that expected from 
a classical treatment. For antiferromagnetic (AFM) ex- 
change coupHngs J < 0, however, for ksT <C |J|, the 
low-T results are very non-classical, even for S — 5/2. 
M(B) for the AFM spin-5 dimer exhibits 25 discrete 
steps, reminiscent of the transverse conductivity in the 



integer quantum Hall effect. In addition, Cy(T, B) ex- 
hibits 25 doublet peaks centered about the corresponding 
magnetization step fields, the splitting of which is pro- 
portional to T. 

At arbitrary T, B, we then evaluate the real part of the 
time autocorrelation function Cii(t), and focus upon its 
Fourier transform Cii(w), which is applicable to inelastic 
neutron scattering experiments. For B = 0, Cii(w) for 
the spin-5 dimer has 45 -I- 1 equally-spaced modes with 
frequencies tj^ „ — n\J\, where n = —25, —25+1, ...,25— 
1,25. Depending upon T and the sign of J, the rela- 
tive importance of these modes varies significantly. For 
B 7^ 0, each of these Ug ^ modes is split into 3 modes. We 
study the T and B dependence of the most important of 
these modes for 5=1/2 and 5/2, for both FM and AFM 
exchange couplings. For the FM case, only a few modes 
are important at low T, and their relative strength is 
nearly independent of \B\. For the AFM case, however, 
the situation is more complicated, as many modes can 
be important at rather low T, and their relative impor- 
tance shifts with \B\. For comparison, we also present the 
analogous results for classical spins. Finally, at arbitrary 
T, B, we evaluate the magnetic dynamical structure fac- 
tor 5(q, u) measurable with inelastic neutron scattering, 
and identify a method by which it can measure Cii{lu). 



THERMODYNAMIC PROPERTIES 
Partition Function 

Here we derive the partition function Z for the quan- 
tum dimer for B 7^ , with spin i = 1 , 2 represented 
by the operator Si, and set H = 1 for convenience. The 
Hamiltonian is 
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FIG. 1: Plots of the static magnetization M normalized to 
2SgfiB versus B = giiBB/\J\ for the AFM S = 1/2 and S = 
5/2 dimers at a = J/{2kBT) = —0.1, —1, —10, as indicated. 



where 7 = (^/i^ is the gyromagnetic ratio. Letting the 
total spin operator s = Si + S2, 7i is rewritten as 



n = -J(s2-S2-S2)/2-7Bs3. 



(2) 



The dimer quantum states are then indexed by the quan- 
tum numbers s and m, where 



s^|sm) = s{s + l)\sm), 
Sz\sm) = mlsm). 



(3) 
(4) 



and since S||sm) — S{S + l)\sm), where S{S + 1) is a 
constant for all measurable quantities, we drop the terms 
proportional to and S| in Ti. for convenience. 
_Letting /3 = l/iksT), a = pj/2, b = p-fB, and 
B = gfiBB/\J\ = 5/(2|a|), the partition function Z = 
Trexp(— /37i) for the Heisenberg dimer is simply 
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Z = J2 J2 e""("+i)+'"" =^I?,(a,6), (5) 
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where 



■Ds{a,b) = 



,+i) Sinh[(2,s + l)b/2] 
sinh(6/2) 



Magnetization 



(6) 



The magnetization M = /3 ^"^elln.^) for the spin-S* 
dimer in a field is then easily found to be 
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M(a,5) = I-Y^V,ia,b)Bsib) 

s=0 

^ lB{Bs{b)), 

where 



(7) 
(8) 



Bsib) = (s + 1/2) coth[(s + 1/2)6] --coth(&/2) 

(9) 

is the standard Brillouin function for a spin s paramag- 
net in a magnetic field. We note that (. . .) represents a 
thermodynamic average over the dimer quantum s val- 
ues, evaluated using the weighting function Vg^a, b). 

In Fig. 1, we plotted M/{2SgfiB) as a function of 
B for the AFM = 5/2 and S* = 1/2 dimers at a = 
-0.1, -1, -10, respectively. At high fc^T/l J| (|a| << 1), 
the dotted curves for a = —0.1, 5 = 1/2, 5/2 are smooth 
functions of B, not too different from the analogous re- 
sults obtained for classical spins. At a = — 1, the dashed 
curves remain smooth, but are shifted over to smaller 
field values, with hints of a kink just below the satura- 
tion magnetization value. Most interesting are the low-T 
effects. At a = — 10, M has 25* steps at integral val- 
ues of i3, as shown by the sohd curves for S — 1/2, 5/2 
in Fig. 1. Since these steps are thermodynamic, they 
are reversible. However, a measurement at fixed B near 
to a step value can lead to interesting, non-monotonic 
behavior of M(T) at low T. 

These thermodynamic steps in M are a consequence 
of quantum level crossing due to the strong B.^ The 
energy Egm of the state \sm) is given by 



-Js(s + 1) 



mB. 



(10) 



where B = ^B and B = B/\J\. For the AFM case, 
J — — |J|, the lowest energy state for each s value is 
Esin,niin — Egs- Thc difference in energy between the 
lowest energy state with quantum number s and that 
with the next highest quantum number, s -I- 1, is then 

A^,, = |J|(s-f 1-|S|), (11) 

which vanishes at \B\ = s + 1. For S = 1/2, there will 
only be one step, as this crossing can only occur between 
states corresponding to the s = and s — 1 quantum 
numbers. Similarly, for S — 5/2, there will be 5 level 
crossings, corresponding to s = 0, . . . , 4. 

For the FM case, the situation is rather boring by com- 
parison. We found that for both S = 1/2,5/2, M{a,b)/j 
is closely approximated by B2s{b)- For a — 0.1, this ap- 
proximation is accurate to a few percent, but for a = 10, 
the corresponding curves for M and B2S are indistin- 
guishable. This is because the s = 2S term in both the 
numerator and the denominator of {Bg) is dominant for 
J > 0. 



Specific Heat 



The specific heat Cy is given by 
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FIG. 2: Plots of Cv/kB versus B at a = -0.1, -1, -3, -10, 
as indicated, for the S = 1/2 AFM dimer. 
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FIG. 3: Plots of Cv/kB versus B at a = -0.1, -1, -3, -10, 
as indicated, for the S = 5/2 AFM dimer. 



leading to 

Cv/kB = 



(s{s + 1)) - coth{b/2){B,ib)) - (Bsib))' 



+2ab 



{sis + l)BM)-{BAb)){s{s + l)) 



(13) 



In numerical evaluations of Cy , we checked that Boltz- 
mann's law, dTCy/T = 2fcs ln(25 + 1), is satisfied. 
As for the FM Af, Cy for the FM spin-5 dimer is rather 
boring. It is closely approximated by the s = 2S terms 
in each of the sums. In this approximation, only the first 
term proportional to b^ survives, so Cy reduces to the 
standard Schottky result for a spin 25 paramagnet, 



Cv/kB C2s(fo), 



(14) 



= 6^(25(25+ 1) + icoth2(6/2) 



-(25* + 1/2)2 coth2[(25 + 1/2)6] ) .(15) 



Numerically, for S = 5/2, the FM Cv is nearly indistin- 
guishable from C2s{b) for a > 1. We note that C2s(&) 
grows as 25(25 + 1)6^/3 for 6 < 1, decays as b"^ exp{-b) 
for 6^1, and has a maximum for b ~ 1. For 5 = 5/2, 
the maximum occurs at 6 « 0.78. 

On the other hand, Cy for the AFM spin-5 dimer 
is much more interesting. In Fig. 2 and 3, we plotted 
Cy/ks versus B at a — -0.1,-1,-3,-10, for dimers 
of spin 5 — 1/2 and 5/2, respectively. The curves for 
a — —0.1 have broad maxima at fields too high to ap- 
pear in these figures. But, as T is lowered to a — —I, 
Fig. 2 illustrates that this broad maximum develops into 
two peaks centered at B w 0,2, respectively. Then, as 
T is lowered further, Cv for 5 = 1/2 becomes two well- 
defined peaks symmetrically centered about B = I, with 
a splitting between them proportional to T. 

In Fig. 3, the analogous AFM Cv for 5 = 5/2 is shown. 
For a = —1, instead of two broad peaks, as in Fig. 2, for 
5 = 5/2 there is an irregular pattern resulting from many 
accessible energy levels. The situation becomes clearer 
at a = —3, with six rather symmetric peaks roughly 
centered at half-integral values of B. However, as for 
5 = 1/2, the low-T limiting behavior is not reached until 
a — —10. For 5 — 5/2, Fig. 3 shows that Cy consists of 
five double peaks symmetrically centered about B = s+1 
for s = 0, . . . , 4, and nearly vanishes at those points. 

We note that this multiplicity of double peaks at low T 
is also a consequence of quantum level crossing. At low 
T, the energy difference between the two lowest states is 
given by Eq. (pi]). At the exact level crossing, these levels 
are degenerate, leading to an exponentially small value of 
Cy, [oc exp(2a)] Just away from these points, excitations 
from the ground state to the first excited state can occur 
for fcsT > I J||s + 1 - \B\\. This should lead to a peak on 
each side of the level crossing points, with a splitting AB 
of the double peaks approximately equal to l/|a|. From 
the data in Figs. 2 and 3 for |a| = 3, 10, the splitting is 
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FIG. 4: Plots of Cv/kB versus kBT/\J\ at B = 
0,0.4, 0.8, 1.0, 1.2, as indicated, for the S = 5/2 AFM dimer. 
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« 2.4/|a|, which is cx T, but 2.4 times as large as in the 
above crude estimate. 

We remark that the 25" double peak structures in Cv 
shown in Figs. 2 and 3 comprise a new quantum effect. 
Although it may be difhcult to increase B to 40 T while 
holding T fixed at ^ 1 K, it might be easier to lower T in 
a fixed, strong B. From Figs. 2 and 3, for B fixed at 0.8 
or 1.2, non-monotonic Cy(T) behavior upon lowered T 
is predicted. However, at S = 1, the behavior at low T 
is very different, as Cv at a = —10 is vanishingly small. 
This striking sensitivity to the precise value of B is pic- 
tured in Fig. 4 for 5 = 5/2. In Fig. 4, for fcsT/l J| > 1, 
the T, B dependencies of Cy are monotonic. Aside from 
the rather ordinary peak for B — 0, the unusual behavior 
is illustrated by comparing the curves for B = 0.8, 1.0, 1.2 
for fcsT/l J| < 0.5 The curves for B = 0.8, 1.2 both have 
peaks at low T, but the curve for the intermediate value 
B = 1.0 does not. This highly sensitive dependence of 
Cv upon T, _B is a new quantum level crossing effect. 



SPIN DYNAMICS 

We now evaluate the time autocorrelation function for 
both quantum and classical spin dimcrs for B 7^ 0. 



Quantum Spin Dynamics 

The time evolution of quantum spins is given by the 
commutator of the spin operator with the Hamiltonian, 

iS,{t) = [S„H], (16) 
S^ = exp(im)S40)exp(-im). (17) 

It is easiest in the quantum case just to keep the time 
dependence in this form, and then to let Ti. operate on 
the eigenstates in the matrix elements of the correlation 
function. The total spin operators and Sz are inde- 
pendent of t. 



Classical Spin Dynamics 

The two classical spins Si (t) each precess according to 
classical Heisenberg dynamics, 

S, = JS, X S + 7S^ X B, (18) 
where |Si(t)| — 1, for i — 1, 2, leading to 



S = 7S X B. 



(19) 



Hence, the total spin S precesses about B, and the Si 
precess about both B and S. The solution to these equa- 
tions is given in the Appendix. 



TIME CORRELATION FUNCTIONS 

Quantum Spins 

In this section, we evaluate the time correlation func- 
tions for quantum spins with general spin value S. Here 
we calculate the time autocorrelation function for the 
quantum dimer. For easy comparison with the classical 
results, we normalize the correlation functions such that 
Re Cii(O) = l.[|9| Then Cii(t) = (Si(i) • Si(0)), or 
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y.e-^^\sm)/D, 

s,m s' ,m' a—1 

x\{sm\Sia\s'm') /D, 

-f3Es„i = as{s + 1) + bm, 
D = S{S+1)Z, 



(20) 



(21) 
(22) 
(23) 



where Z is given by Eq. (||). In Eq. (pT|), the expectation 
values are related to Clebsch- Gordon coefficients. 

For quantum spins, the near- neighbor correlation func- 
tion and the autocorrelation function are related by the 
sum rule, 

{sl) + {s'-sl)cos{Bt) - i(s(t).s(0) + s(0) •«(*)) 

= 2S'(S'+1) 

xRc[Ci2(t)+Cn(i)], (24) 

where Ci2(i) is normalized as is Cii{t). The thermal aver- 
ages (s^) and (s^) can be written in the notation of Eqs. 



(|6|)-(|9D as {s(s + 1)) and (s(s + 1)) - coth(6/2)(B46)), 
respectively. 



Classical Spins 

It is interesting to compare the quantum time corre- 
lation functions with those obtained from classical spin 
dynamics. For classical spin dimers, one solves the equa- 
tions of motion of the spins directly, taking the length of 
each spin to be unity. We first solve Eq. (|l9|) for S(t), and 
then Eq. ( |l8| ) to obtain the Si(t). We then evaluate the 
classical autocorrelation function Cii(i) — (Si(i) • Si(0)) 
by averaging over the length S* (0 < 5 < 2) of the to- 
tal spin, the angle 9 between S and B, and the angle 
describing the initial relative configuration of Si and 
S2, and requiring Cii(O) = 1. Since the procedure is 
analogous to that used for the four-spin ring,|T^ the re- 
sults are given in the Appendix. Here we only remark 
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that the classical correlation functions also must satisfy 
a sum rule, 

l[l + 2cos{Bt)]{S^) = Cn{t)+Ci2{t). (25) 



FOURIER TRANSFORMS 



Quantum Spins 

The Fourier transform Cii{lu) of the real part of the 
autocorrelation function Cii(i) for the quantum dimer of 
spin- 5 spins is given by 

6S+1 

Cn{uj) = J2 fsASi'^ - ^Si) + S{u; + ujsi)]. (26) 

i=0 

The discrete mode frequencies lusa and their ampli- 
tudes fs,i are given for S — 1/2 and S = 5/2 in Tables 
II and III in the Appendix. 



Classical Spins 

In order to compare the quantum and classical re- 
sults, it is useful to evaluate the Fourier transform 
Cii{uj) = J^^dtcxp{iuit)Cii{t) of the classical autocor- 
relation function. Letting w = co/\J\ and B = B/\J\, we 
find for positive w that 

Cn(w) = S{Q)Coo + S{u-B)Coi+SCn{i:')- (27) 

Exact expressions for the discrete amplitudes Coo ^md 
Coi, and for the continuous part 6Cii{lu) are given in the 
Appendix. Then, in order to check the numerical eval- 
uation of the above quantities, we employ the frequency 
sum rule, 1 = /q (iti)Cii(a;)/7r. 

We note that Coo also appears in the expression for 
Cii{uj) with a; < 0. The coefficients Coo and Coi of the 
delta functions at w = and uj = B, respectively, are 
functions of both a and B. These quantities arise from 
the long-time rms limit of Cii(f) and for the long-time os- 
cillatory behavior at the driving frequency, respectively. 
This second delta function occurs at the resonant fre- 
quency in an EPR experiment. 



NUMERICAL RESULTS 
Autocorrelation Function Spectra 

In Fig. 5, we plot the amplitudes /5/2,i of the delta- 
function modes of the zero-field quantum 5 = 5/2 AFM 

Cii{lu) for positive ti-'/|J|, at the temperatures T = 
kBT/[\J\S{S + 1)] = 0.1,1,00. The amplitudes of the 
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FIG. 5: Plots of the AFM quantum S = 5/2 modes of 

Cii versus uj/\J\ at B = at T = kBT/[\J\S{S + 1)] = 
0.1,1,00, and for the combined discrete classical modes at 
u! = 0. The curves represent 5Cii{aj)/[S{S + 1)]^''^ versus 
a;[5(5+ 1)11/71 J|. 



modes at these three T values are indicated by the sym- 
bols -|-, X, *. For comparison, we also show the results 
of the zero-field classical calculation Cii(a;). For B = 0, 
both of the delta-function amplitudes Coo and Cqi appear 
at w = 0, and their combined weight at these T values is 
indicated by the circle, triangle, and square, respectively. 
In addition, the continuous part SCii{uj) is plotted both 
by scaling its amplitude by 1/[S{S + 1)]^^^ and the fre- 
quency by [6'(S' -|- 1)]^/^. This combined scaling allows 
us to compare with the quantum results, while preserv- 
ing the area under the curves. For S = 5/2, this scaling 
changes the maximum u) = u)/\J\ from 2 to \/35 w 5.91. 
We note that as T ^ oo, the quantum and classical delta 
functions at a; = agree exactly, both having the weight 
7r/2. Also as T — > oo, the classical ^Cii(a;) forms an enve- 
lope for the amplitudes of the quantum delta functions, 
except for a slight deviation at the larger frequencies. At 
T = 1, the agreement between the quantum and classical 
results is also pretty good, although there are deviations 
at nearly every quantum mode value. These deviations 
are more pronounced at T = 0.1, but the overall agree- 
ment is still rather good. 

In Fig. 6, we show the corresponding curves for 
the 5 = 5/2 AFM dimer at the field B = 0.4. For 
-B 7^ 0, each of the quantum zero-field modes splits 
into three modes, and therefore to compare the contin- 
uous part (5Cii(a;) with the split quantum modes, we 
divide the scaled continuous part of Cn by 3, plotting 
5Cii/3[S{S+l)]^/^ versus w[S'(5-H)]i/V|J|. As in Fig. 
5, the classical delta functions modes at cD = 0, B are not 
scaled, but in this case, they are distinct. For this rather 
strong field value, B = 0.4, the agreement between the 
classical and quantum results is very good as T — > oo, 
and pretty good at the intermediate T = 1. However, at 
T = 0.1, these results can differ by more than a factor 
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classical quantum T/[|J|S(S+1)] 
(scaled) values 
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FIG. 6: Plots ofjhe AFM quantum 5 = 5/2 modes of Cn 
versus a;/|J| at _B = 0.4 at T = 0.1,l,oo, and for the dis- 
crete classical modes at cj/|J| = 0,B. The curves represent 
SCii/3[S{S + versus lo[S{S + 1)Y^'^/\J\. 
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FIG. 7: Plots of_the FM quantum 5 = 5/2 modes of Cn 
versus i^/\J\ at B = at T = 0.1, 1, oo, and for the com- 
bined discrete classical modes at a; = 0. The curves represent 



excellent as T ^ oo, pretty good at T = 1, but only 
qualitative at T = 0.1. In this case, the classical delta 
function at = differs significantly from the quantum 
one, and the continuous classical curve cannot be scaled 
in this way at all. As shown for the equivalent neighbor 
model, the FM classical curves obey a different low-T 
scaling relation, as the low T peak center approaches NJ 
linearly in T from below as T ^ 0.|l| In Fig. 8, the 
development of the continuous classical curve into three 
low T modes is evident, but the fit to the quantum case 
at T = 0.1 is not good at all. Not only should the low-T 
classical curves require a different scaling relation, similar 
to those for B = 0, the order of the relative heights of the 
three largest modes at the largest frequencies pictured 
for the classical curves is inverted relative to that for the 
quantum case. 

We now consider the case S — 1/2. For the B = 
AFM and FM cases, there are two quantum modes for 
< w, at uj/\J\ = 0, 1. As T oo, these modes have 
equal intensity (7r/2), for both FM and AFM cases. For 
the AFM case as T decreases to 0, the amplitudes of the 
modes at uj/\J\ — 0(1) decrease (increase) continuously 
to (tt), respectively. For the FM case, however, these 
mode amplitude values increase (decrease) with decreas- 
ing T to 27r/3 (tt/S), respectively, at T = 0. Somewhat 
surprisingly, for the AFM case, the amplitude of the com- 
bined classical delta functions at w = closely tracks 
that of the quantum oj = mode. For the FM case, how- 
ever, this tracking is not so good. However, the best that 
can be said for the continuous part is that as T ^ oo, 
the second quantum mode nearly falls upon the classi- 
cal curve scaled as in Figs. 5 and 7, for both AFM and 
FM cases. At T = kBT/[\J\S{S +1)] = 1,0.1, this agree- 
ment becomes increasingly much worse, respectively. For 
the FM case, the classical treatment fails miserably as 
T ^ 0, as the amplitude of the classical delta function 



<5Cii(a;)/[S(5' + 1)]^/^ versus LJ [5(5'- 



11/2 



/\J\- 



of two. We note that the classical curves are piecewise 
continuous functions of w, with several discontinuities in 
slope evident. These discontinuities in slope generally ap- 
pear at the the frequencies corresponding to rather large 
jumps in the quantum mode amplitude values. At T val- 
ues lower than those pictured, the classical curves develop 
into three distinct modes with finite widths. 

In Figs. 7 and 8, we show the results for the FM S = 
5/2 case a.t B = 0,0.1, respectively. Choosing B = 0.1 
in Fig. 8 allows us to display all of the results clearly 
on the same figure. Except for the different T values in 
Figs. 8 and 6, Figs. 7 and 8 correspond precisely to 
Figs. 5 and 6, with the same scaling of the continuous 
part of the Fourier transform of the classical correlation 
function, and the same symbols and linestyles. For the 
FM case at B = 0,0.1, Figs. 7 and 8 show that the 
agreement between the classical and quantum modes is 
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FIG. 8; Plots of_the FM quantum 5 = 5/2 modes of Cn 
versus ijj/\J\ at B = 0.1 at T = 0.1, l,oo, and for the dis- 
crete classical modes at li'/IJI = 0,B. The curves represent 
5Cn/3[5(5 + l)]i/2 ygjg^g ^[5.(5. ^ 1)]^^V\J\- 
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FIG. 9: Plots ofjhe AFM quantum S = 1/2 modes of Cii 

versus a;/|Jl at B = 0.4 at T = 0.1, 1, oo, and for the dis- 
crete classical modes at u)/\J\ = 0,B. The curves represent 
5Cii/3[S{S + l)]i/2 versus u;[S{S + 1)]'/V| J|. 
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FIG. 11: Plot^of the FM quantum S =1/2 mode amplitudes 

f 1/2,1 versus B at T = 0.1 The numerical labels correspond 
to the i values in Table I of the Appendix. 



mode at w = approaches tt, and the integrated inten- 
sity of the continuous classical mode peak vanishes as 
T ^ 0. Hence, the classical treatment does not describe 
the 5=1/2 dimer at i? = correctly, and the correct 
quantum treatment leads to just two modes for w > 0. 

For B ^ Q, however, there are five modes for a; > 
with 5 = 1/2, so a classical treatment can approximate 
the quantum behavior a bit better than for _B = 0. In 
Figs. 9 and 10, the AFM and FM 5 = 1/2 results for 
B = 0.4 at T = 0.1, 1,00 are shown. In both cases, the 
two quantum modes for w/IJ] =0,5 = 0.4 are com- 
pared with the two classical delta functions, and the re- 
maining three quantum modes for < a; are compared 
with the continuous classical curves. In both cases, as 
T ^ oo, the five quantum modes have the amplitudes 
7r/6, 7r/3, 7r/6, 7r/6, 7r/6, respectively. The classical results 



agree precisely with the first two, and form a qualitative 
envelope similar to the remaining three quantum modes. 
This agreement is qualitatively preserved at T = 1. How- 
ever, as T ^ 0, the classical and quantum cases differ 
dramatically. For the AFM case pictured in Fig. 9, the 
five quantum modes for < a; approach the amplitudes 
0, 0, 7r/3, 7r/3, 7r/3 as T — > 0, respectively, and the con- 
tinuous classical curve develops a strong peak at w = 0, 
which is qualitatively different. The qualitative behav- 
ior of the two discrete classical AFM modes is not too 
bad, however. For the FM case pictured in Fig. 10, 
as T 0, the quantum FM modes approach the ampli- 
tudes 7r/3, 7r/3, 0, 0, 7r/3, respectively. That is, the modes 
at w = 0, 1 J -|- B\ increase, the modes at a; = | J|, | J — B| 
decrease, and the mode dXu) = B has a non-monotonic T 
dependence. The classical treatment preserves these FM 
features only qualitatively, and is inaccurate for T < 0.1. 
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FIG. 10: Plots of the FM qu_antum 5 = 1/2 modes of Cn 
versus a;/|J| at _B = 0.4 at T = 0.1, 1, oo, and for the dis- 
crete classical modes at u}/\J\ = 0,B. The curves represent 
5Cii/3[S{S + 1)]^/^ versus u;[S{S + 1)]'/V| J|. 



Field dependencies of the quantum modes 

In Fig. 11, we plotted the amplitudes /i/2,i of the five 
FM quantum modes for = 1/2 at T = 0.1 versus B. 
The modes 0^1/2, i for i = 0,4 correspond to a; = 0, \B\, 
\J\, \ J — B\, and |J-|-i3|, respectively. Note that at -B = 
0, the u)i/2^i for i = 0,1 are degenerate but unequal in 
intensity, and the L0i/2.i for z = 2, 3, 4 arc both degenerate 
and equal in intensity. At each field, the sum of the 
five intensities sums to n. From this figure, it is evident 
that the crossover from the low-field regime to the high 
field regime occurs at the rather low field value, B « 
0.1. The high field regime is clearly consistent with the 
mode amplitudes at B = 0.4 pictured in Fig. 10, with 
the amplitudes of the modes at i^/\J\ — 0,0.4,0.6, 1, 1.4 
approximately equal to 7r/3, 7r/3, 0, 0, 7r/3, respectively. 

In Fig. 12, we plotted the amplitudes fi/2.i of the five 
AFM quantum modes for S* = 1/2 at the T = 0.1. From 
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Fig. 12, the five modes arc difficult to discern clearly, 
due to the strong degeneracies involved. Clearly, the 
crossover from the low-field to the high-field regime oc- 
curs at i? = 1. In the low- field regime, the modes at 
oj = 0,B have nearly zero amplitudes, and the modes 
with u) = I J|, I J ± B\ have nearly equal amplitude n/3. 
This is the situation pictured for B = 0.4 in Fig. 9. How- 
ever, Fig. 12 indicates that dramatic changes in the mode 
amplitudes at low T should take place as B is increased 
from w 0.8 to « 1.2. Over this field range, the amplitudes 
of the modes at w = | J|, | J + i3| decrease from nearly tt/3 
to nearly 0, and the amplitudes of the modes at a; = 0, B 
increase from nearly zero to nearly 7r/3. Meanwhile, the 
amplitude of the remaining mode at w = \J\ remains 
constant at w/S. Although not pictured for brevity, at 
T = 1, all five AFM modes for S* = 1/2 are clearly evi- 
dent, and the crossover from the weak field regime to the 
strong field regime occurs at B ~ 1.5. 

The field dependencies of the most important modes 
for the FM and AFM 5 = 5/2 dimers are shown at the 
rather low temperature T = 0.1 in Figs. 13 and 14, 
respectively. In Fig. 13, the FM modes shown are for 
LO = 0, B, \5J±B\, and5|J|. The crossover from the weak 
to the strong field limits occurs at about B 0.5. In the 
weak field limit, this corresponds to the dominant modes 
at B = 0.1 shown in Fig. 8, for which the ranking of the 
five strongest low-T modes is at w/| J| = 0.1, 0, 5.1, 5, 4.9, 
respectively. At high fields, B > 2, there are only three 
important modes, at w = 0, B, and 5 J -|- B, which have 
amplitudes that approach 57r/7, 7r/7, 7r/7, at B — > oo 
respectively. 

For the AFM case with S = 5/2, the behavior of the 
relative mode amplitudes at T = 0.1 is more complicated, 
as shown in Fig. 14. There we plotted the field depen- 
dencies of the eight most important modes, which occur 
at w = 0, |2J±B|, 2\J\, \B\, |J±S|, |J|, respectively. 
Prom Fig. 14, it is evident that there are essentially three 
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FIG. 12: Plots of th_e AFM quantum 5=1/2 mode ampli- 
tudes /i/2,i versus B at T = 0.1. The mode i values corre- 
spond to Table I in the Appendix. 
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FIG. 13: Plots of the dominant_FM quantum 5 = 5/2 mode 

amplitudes fr,/2,i versus _B at T = 0.1. The mode i values 
correspond to table II in the Appendix. 



field regimes. The low field regime occurs for B < 0.5, 

the intermediate regime is for 1.5 < B < 2.2, and the 
high field regime for _B > 3, approximately. An exam- 
ple of the low-field regime at the same temperature was 
shown for B = 0.4 in Fig. 6. In Figs. 6 and 14, the 
low-field rank-ordering of the six most important modes 
are those at w = \B + J\, |J|, \B-J\, \B + 2J\, 2\J\, and 
\B — 2J\, where J < 0. From the intermediate regime 
pictured in Fig. 14, however, the two most important 
modes are those at ui = \B + 2J\ and \B + J|, with the 
six other modes in the list contributing lesser, but com- 
parable amounts. Finally, in the high-field regime, the 
two dominant modes are at w = and \B\, respectively. 
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FIG. 14: Plots of the dominant AFM quantum S = 5/2 mode 
amplitudes /5/2,i versus B at T = 0.1. The mode i values 
correspond to table II in the Appendix. 
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DISCUSSION 

We remark that there have been some experiments on 
one or more of the Fe2 single molecule magnet dimers. Le 
Gall et al. measured the four Fe2 dimers they made, and 
found that the zero-field magnetic susceptibility could be 
fit with the Heisenberg model, with exchange constants 
ranging from 15 to 19 cm~^ (22 to 27 K, which corre- 
sponds to B = 1 at S « 16 - 22 T). A sHghtly different 
fit was made by Lascialfari et al, obtaining |J| w 22^ 
for the Fe2 dimer, [Fe(OMe)(dpm)2]2- lIHI Those au- 
thors also refer to EPR measurements that give rise to 
a zero-field splitting anisotropy of about 7 K in the first 
excited state, which would complicate the analysis per- 
formed here. Such a zero-field splitting would still allow 
quantum level splitting, but the magnetization steps and 
the specific heat double peaks would not all be regularly 
spaced. We remark that these level crossing effects are 
purely quantum in nature, as the analogous classical cal- 
culations do not give rise to such effects, and hence are 
only approximate for \a\ ^ 1. Nevertheless, the quantum 
effects predicted for the magnetization and specific heat 
ought to be observable with presently available facilities, 
as one should be able to see one or two magnetization 
steps, and one or two specific heat double peaks. 

In addition, inelastic neutron scattering would see 
slightly different results from those presented here. In 
this situation, the classical envelope curves might be 
quite useful, as they can serve as a guide to the behavior 
that might be expected with specific zero-field splitting 
values. More important, inelastic neutron scattering can 
be used to probe the details of the magnetic interactions 
within an individual dimer. Provided that a single crys- 
tal of sufficient size for such studies can be obtained, 
one would perform the experiments at the wave vectors 
specific to the crystal structure. More precisely, the un- 
polarized inelastic neutron magnetic dynamical structure 
factor S{q, to) for a single crystal of equal-spin dimers is 
given by I I H 



.13=1 

y ^^e-*(Ql(q,i)Q^(q,0)), (28) 



Q,/3 = l 

X 



where the qa = sin 9q cos (j)q , sin 9q sin (j)q , cos dq for a — 
1,2,3, respectively, and 9q, (j)q are the angles the scat- 
tering wave vector q makes with the spin quantization 
coordinates, and (. . .) represents a thermal average, as in 
Eq. (pO|). The dimer structure operator 



Qa(q,i) = /(q)[e*''-''5ia(t) + e''^i"^52a(i)], (29) 

where 2d is the displacement vector between the dimer 
spins, and /(q) is the atomic magnetic form factor. 
Pi The low T, B = transitions in 5 = 1/2 Cn^+ and 



Yb^"*" and in S* = 3/2 Cr^+ dimer powders were treated 
previously. H, ^ |l^ For B = 0, the quantization axis is 
arbitrary, so 9q = n/2 was chosen.]^ For B 7^ 0, the 
quantization axis is parallel to B, so 9q and (j)q are the 
angles q makes with the coordinates z| |B and x x y = z. 
For scattering wavevectors qc directed along the special 
angle 

9^ = sin-i (2/3)1/2, (30) 

the components of each Cij{uj) are sampled equally, and 
2 



37r 
2 



/2(q,)[Cii(w) -|-Ci2(t^)cos(2q, • d)] 

(31) 



= |/'(qc) Cii(^)[l-cos(2qe-d)]/7r- 



cos(2qc 



S{S+1) V"''--'"' ' 2 
x[S{lj - B) + 5{u; + B)] 



(32) 



where we have employed the Fourier transform of the sum 
rule in Eq. (|2^). For powder samples, one can still use 
the special angle technique with a field to obtain Cii(a;), 
but since the direction of d is random, one obtains 



cos(2qc • d) 



sm(2qcd) 
2qcd 



(33) 



where • • • is a spatial average. p2[ . 

For the general case of 9q ^ 9q, however, S{c[,uj) can- 
not be written simply in terms of the C,;j(w). There are 
four factors hi for i — 0, . . . , 3, listed in the Appendix, 
that depend upon q • d and sin6'„. We then find 



27r 



fs.ihs,i 



(34) 



where hs^i is the appropriate hi for the modes iwsi, as 
indicated for S — 1/2 and 5/2 in the Appendix. The 
factors /lo and hi are cx cos2(q • d), and correspond re- 
spectively to the w = and ±5 modes. The factors /i2 
and /13 are oc sin2(q • d), and correspond respectively to 
the modes at ±nJ and ±|nJ± B|, for n— 1, . . . , 25 + 1. 
Since /iq and /12 are also (x sin^^g, whereas hi and h^ 
are cx (1 — ^ sin^ 9q)^ the experimenter can fine tune the 
single crystal data by rotating B and q relative to d. 

Neutron powder data on the deuterated 5=1/2 dimer 
single molecule magnet VODP04-iD20 (V2) were taken, 
resulting in a fit to the AFM Heisenberg model of | J| = 
7.81(4) meV, Q close to the value 7.6 meV found 
the susceptibility fit 



1 20 



m 
at 



This corresponds to i? = 1 
S « 66 T, which is too large for thermodynamic studies. 
However, inelastic neutron scattering at 0.1-0.25 ought 
to be possible for this material. Single crystal data could 



10 



be particularly interesting. For the Fe2 dimers, inelastic 
neutron scattering in a field of i? ^ 0.4 — 0.5 should be 
possible, which would not see any level crossing effects, 
but could prove interesting, as indicated in Fig. 6. 



Fi{x) = F^{x) = [cosh(a;) — sinh(x)/a;]/a; and ^2(0;) = 
Fq{x) = Fo{x) — 2Fi{x)/x, respectively. We note that 
Fi{x)/Fo{x) — L{x) = coth(a;) — 1/a; is the Langevin 
function. Wc find 
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APPENDIX 
Classical Time Correlation Function 

From the classical equations of motion, 

S(i) = BS"!! + S'i[icos(i?<) - ysin(Si)], (35) 

where B = 7_B, B is a unit vector parallel to B, and x, 
and y are orthogonal unit vectors satisfying x x y = B. 
Since 6*11 = S • B, we have S\\ — Scosd and S± = S'sin6', 



and hence = Sf, + S^. Si{t) is then found to be 



^1,11 W = ^1,11 +5i,±cos(J^t-0o), 



(36) 



Si±it) 



Si^\\S±_ Si^j_S± 



'exp[^i{JSt - (I)q)] exp[±i(JS't - (/)o)] 



S 



(37) 



where Si± = Six i iSiy^ and (/)o fixes the initial relative 
configuration of the two spins. Since S2 = S— Si, we have 
Si,ii = 82,11 = ^11/2 and Si,i_ - -S2.± = iS^/S)[l - 
5'2^4]i/2^ After averaging over ^o, 



Cii(t) 



2^(1 



+Sl cos{Bt) 



-574) cos{SJt) 
(1 - 574) 



52 



X coa[{JS + B)t] + {S - S^i)^ 
xcos[{JS - B)t]]). (38) 



Replacing S with s for elegance, these classical averages 
are evaluated from 



Z 



Z-^ / sds 
Jo Jo 

sds 

10 JO 



,QS COS £ 



sm 



as -i-bs cos 9 



(39) 



(40) 



The integrals over 6 can then be written in terms of 
Fq{x) = sinh(a;)/a; and its first and second derivatives. 



- 8Z 



sdse° 



2s^Gi{bs,Bt) 



+ {4- s-'){cos{sJt)G2{bs,Bt) 



+G3(6s,Si)sin(sJt)jJ, (41) 

Gi{x,y) = F2{x) + [Foix)-F2ix)]cos{y), (42) 
G2{x,y) = Fo{x)-F2{x) + [Fo{x)+F2{x)]cos{y), 

(43) 

G3{x,y) = 2Fi{x)sm{y). 



As r ^ 00. we set t* 



\J\t, and obtain. 



[lJ^2cosCBi)] 



fit* 



1 



[1 + 2cos(2r)] 
3sin(2i*) 3[1 - cos(2i*)] 



t* 



2t* 



(44) 



(45) 



(46) 



We note that f{t*) was obtained previously for the 
zero field case. 1 13 We also have limT^oo Ci2(t) = 1 — 
limT^^oo Cii(t). These forms clearly satisfy the require- 
ment limy^cx) Cii(O) = 1. 



Quantum Frequency Spectrum 

The Fourier transform Cii(ti;) of the real part of the 
autocorrelation function Cii(i) for the quantum dimer of 
spin-S spins is given by 

6S+1 

Cll(cj) = ^ fs,i[S{^-^Si)+5{LO + LOS^)], (47) 

2 = 

fs,^ = ^as7C^(&,a). (48) 
For both 5 = 1/2 and 5 5/2, C-^(5, a) is given by 

C^ib.a) = 45(5+l)e2S[''~(25+i)"lz. (49) 

We note that D = S{S+l)Z and that C^(0,0) = 45(5 + 
1)(25+1)2. 

The factors hi that weight the modes in 5(q, uj) given 
by Eq. ( ^ can be derived from Eqs. (|8|) and (|2|). In 
Eq. (28|), the off-diagonal terms in 5(q, w) with a, f3 = 
1, 2 sum to zero, and the remaining off-diagonal terms all 
vanish. Hence, we only require the matrix elements 



M 



s m 
sm,a 



= (l-ql) (5m|Qi(q,0)|.W) //7q).(50) 
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i 


As 


Am 


hi 











sin^ 6q cos'^(q ■ d) 


1 





±1 


(1 - isin^ 6l,)cos2(q-d) 


2 


±1 





sin^ 6q sin^(q ■ d) 


3 


±1 


±1 


(1 - isin^ 6'Jsin2(q. d) 



TABLE I: Factors hi that appear in S{q,uj), Eq. (31), and 
their associated transition quantum number changes. 



i 


^1/2, i 


h\/2,i 


«l/2.i 








ho 


l + e^' 


1 


\B\ 


hi 


[l + e'f 


2 


\J\ 


h2 




3 


\J-B\ 


hs 


e''(e' + e-2") 


4 


\J + B\ 


hs 


l + e''-^- 



TABLE IL Frequency us.i and weighting factor hs,i spectra 
of the 5 = 1/2 Heisenberg dimer in a magnetic field. 



From Eq. ^ and s±\sm) — Af^\s, m± 1), where Af^^ — 
[s{s + 1) - m(m ± 1)]^/^, we write (|| 



{sm\Snz\s'm' 
{sm\Sn±\s'm' 



x{BsmSs',s + l + CsmSs',s-l)], (51) 



= (5,„',,„±i[AJ,„(5s'.s/2 + (-1)" ^ 

x{Bt„,Ss'^s+l+CtmSs',s-l)], (52) 



for n = 1,2, respectively. We note that 
s± = 5i± + S'2±. We find 



^1 



^2 



\h,{At^YSs>^s/2 + 2h, 



:(0 



Ss',s-1 



,(53) 



where the /i^ and the changes As = s' — s and Am = m' — 
771 in the matrix elements for which they occur are listed 
in Table I. For arbitrary S, each of the hsi is equal to one 
of the four h, hsted in Table I. For 5=1/2 and 5/2, the 
appropriate choices of the hsi are listed in Tables II and 
III, respectively. We remark that by setting /(q) = 1 
and each of the hi = 1, ^(q, w) Cii{lo)/ {2i:). 



Quantum Frequencies for S = 1/2 

For S = 1/2, we have 

C^/^{b,a) = 3(l + e^ + e2'' + e^-2"), (54) 

The mode frequencies Wi/2,i, their relative amplitudes 
ai/2,i, and the factors /ii/2,i appearing in S{q,uj) are 
given in Table II. 



i 


'^5/2, i 


h5/2,i 


0-5/2,1 








ho 


2g5{6-6Q) y^2S ^an(n + l)^ 

X cosh(s6) 


1 




h2 


:|(lVe-2-)e5''-28- 


2 


2\J\ 


h2 


'§e^'~''-{l + e-'-)Y.UAs{h) 
f^e'^'^^^il + e-^nHl-oMb) 

l{l + e'''^-)e^'Y.UMh) 


3 


3|J| 


h2 


4 


4\J\ 


h2 


5 


5\J\ 


h2 


6 


\B\ 


hi 


2 coth(6/2)e5(''-«") Y.^'n^i X 
^ YJl-i'^^'^'^Hsb) 


7 


\B + J\ 


hi 




8 


\B + 2J\ 


h3 


32g3(,-24.(^^gi,-4.)^^(^) 


9 


|B + 3J| 


hi 


||e2'-i»"(l + e''-«-)X3(fe) 


10 


|S + 4J| 


hs 


i|Oef-io^(l + e''-«")X4(6) 


11 


|B + 5J| 


hs 


i(l + e''-i«")X5(6) 


12 


\B~J\ 


hi 


f e5^-2»"(e'' + e-2")Yi(fe) 


13 


\B - 2J\ 


hs 


f2e«''-24"(e'' + e-4")F2(fe) 


14 


|B-3J| 


hi 




15 


|B-4J| 


hi 


140g8i,-10a(g6^g-8«)5^^(^) 


16 


\B - 5J| 


hi 


|e9f(e'' + e-w«)Y5(6). 



TABLE III; Frequency and weighting factor hs.i spectra 
of the S = 5/2 Heisenberg dimer in a magnetic field. 



Quantum Frequencies for 5 = 5/2 

For simplicity, we set 

(2s + l)sinh[(2s + 1)6/2] 



A,{b) 
Xs{b) 



sinh(6/2) 



(55) 



e/'' sinh(s6) 



4sinh^(6/2) 
-ssinh(fe) -4s2 sinh^(V2) 



(56) 



n(6) = X,(-6) = 



4sinh^(6/2) 
ssinh(6) +4s2 sinh^(V2) 



sinh(s6) 



(57) 



We note that Xi{b) = Yi(b) = 1 and that ^^(0) = 
Ys{0) — s(4s^ — l)/3. The mode frequencies ^5/2,1, 
their relative amplitudes 05/2,1, and the h^/2,i factors are 
given in Table III. In evaluating the coefficients a^j2,i at 
a = 6 = 0, it is useful to employ the relations 



2S 



n 

^(2s + l)^ 



5(5 + l)(25+l)2 



(77 + l)(2r7+ l)(277 + 3) 



(58) 



(59) 



s=0 
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Classical Frequency Spectrum 

From Eq. ( p7| ) and letting a) — lo/ J and B = B/J, the 
classical spin Fourier transform Cii{lo) has the following 
discrete and continuous contributions, 



n2{B) = se(2-B) + 2e(B-2), (72) 

n^iB) = 2Be{2-B) + (B + 2)e(B~2). (73) 
In Qi{B), the indicates that the maximum of the mode 
is at (I> = 0, the same position as for Cqo- 



Cii(w) = 

SCii{lu) = 

Coo = 

Cqi — 

Ci{Cj) = 

C2{Co) = 



4 

1 



1=1 



AZ 



e(2-w)2(I'(l -c;)74)e"'^ X 



AZ\J\ 
x[Fo{bCj) ~ F2ibCb)] 



(60) 
(61) 

(62) 
(63) 

(64) 



efremov@m£i£ks^dresden^^^ 



rklemmQmpipks-dresden . mpg. de 



AZ\J\ 



6(2 - S - (I;)(J) + S)[l - (w + Bf/A] x 



xe 



Q(li + B)' 



' (^FoibiC^ + B)] + 



F2[biu + B)] + 2F,[\b\i[b + B)] 



^ -Q{B - Cj)Q{Cj + 2 - B){B - Oj) X 



(65) 



AZ\J\ 

X [1 - (c:; - B) V4]e"(^-^)' (Fo[6(S - Cj)] 
+F2[6(S - w)] + 2Fi[|6|(B - (i)]' 



(66) 



C4(w) 



AZ\J\ 



e(w - 5)6(2 + s - j;)((:;} - s) X 



xe 



a(tD-_B) 



'(fo[6(c:;-s)] 



+F2[6((I'-S)] -2Fi[|6|(w-B)] 



Low Temperature Classical Modes 



(67) 



We now investigate the low T behavior of the various 
contributions Ci{ijj) to the classical Cii(a;). We follow the 
procedure used for the isosceles triangle and equivalent 
neighbor models in zero field. (16[[18| 



The FM modal spectrum as T — > is given by 
n^{B) = \B-2\, 

n2(B) = 2, 

n^iB) = B + 2. 



(68) 
(69) 
(70) 
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We note that Ci leads to 512, C2 and C3 combine to create 
Oi, and C4 leads to ft^. 

As T — > 0, the AFM mode frequencies satisfy 



ni{B) = 09(2 -B) + (5-2)6(5-2) 



(71) 



